UCI-HEP-TR-2010-15 



Constraints on Dark Matter from Colliders 

Jessica Goodman, Masahiro Ibe, Arvind Rajaraman, 
William Shepherd, Tim M.P. Tait and Hai-Bo Yu 

Department of Physics and Astronomy, 
University of California, Irvine, California 92697 
(Dated: April 5, 2011) 

Abstract 

We show that colliders can impose strong constraints on models of dark matter, in particular 
when the dark matter is light. We analyze models where the dark matter is a fermion or scalar 
interacting with quarks and/or gluons through an effective theory containing higher dimensional 
operators which represent heavier states that have been integrated out of the effective field theory. 
We determine bounds from existing Tevatron searches for monojets as well as expected LHC 
reaches for a discovery. We find that colliders can provide information which is complementary or 
in some cases even superior to experiments searching for direct detection of dark matter through 
its scattering with nuclei. In particular, both the Tevatron and the LHC can outperform spin 
dependent searches by an order of magnitude or better over much of parameter space, and if the 
dark matter couples mainly to gluons, the LHC can place bounds superior to any spin independent 
search. 
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I. INTRODUCTION 



While astrophysical observations provide compelling proof for the existence of a non- 
baryonic dark component to the Universe and precise measurements as to its abundance 
lj], they offer no clue as to the mass of dark matter (DM) particles, how they fit into 
the Standard Model (SM) of particle physics, or even whether or not the dark matter has 
interactions beyond gravitational. The most compelling vision of dark matter is a weakly 
interacting massive particle (WIMP), which offers the possibility to understand the relic 
abundance of dark matter as a natural consequence of the thermal history of the Universe 
through the WIMP (less) miracle j^j. The large interactions of WIMPs with SM particles 
may imply detectable rates of WIMP annihilations into SM final states, scattering of WIMPs 
with heavy nuclei, and production of WIMPs in high energy reactions of SM particles at 
colliders. 

WIMPs which produce signals in direct detection experiments must also couple to nu- 
cleons, and therefore can be produced at colliders like the Tevatron and LHC. Low mass 
particles are particularly amenable to searches at colliders, since the fact that a typical colli- 
sion involves quarks and/or gluons carrying only a small fraction of the parent (anti-)proton 
energy implies that cross sections fall dramatically with the mass of produced states. Light 
states can thus be produced with very large rates. In the case of a WIMP, stability on 
the order of the lifetime of the Universe implies that pair production must highly dominate 
over single production, and precludes the WIMP from decaying within the detector volume. 
WIMPs therefore appear as missing energy, and can potentially be observed by searching 
for visible particles recoiling against dark matter particles UQ. This can be used to set 
constraints on the WIMP couplings to the constituents of nuclei, which in turn can be trans- 

n 

lated to constraints on direct detection cross sections. In previous work |8j, this was done 
for the case of Majorana WIMPs 1 ; here we extend this work to also include Dirac fermion 
and scalar (real or complex) WIMPs. 

There is currently particular interest in light dark matter particles. The DAMA ex- 



periment has reported a signal of annual modulation at a high significance level 



3; 



this 



signal is consistent with a dark matter discovery interpretation with a dark matter parti- 



1 A partial set of operators for a Dirac WIMP were bounded using Tevatron data in Ref. [9( and the 
Tevatron limits and LHC reach for operator D8 (see below) were considered in Ref. 
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cle of mass < 10 GeV 
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12]. The CoGeNT collaboration has also reported a signal \U 



which can be explained by a WIMP in this mass range (though there may be some tension 
with unpublished data from 5 towers of CDMS Si detectors [14J , and with recent data from 



XENON10/100 
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161]) . There has hence been much recent interest in models of light dark 



matter (where the DM mass is order a few GeV) 17H26I]. As colliders are most effective 
when producing highly boosted, light WIMPs, the tantalizing hints from DAMA and Co- 
GeNT point toward a region where colliders can have a particular impact on theories of dark 
matter. 

In this article, we consider the situation where the WIMP is the only new particle in the 
energy ranges relevant for current experiments. Given the small energy transfers involved 
in direct detection, this assumption is almost certainly justified. For colliders, the degree to 
which it is justified depends on the details of the WIMP theory. Under this assumption, the 
WIMP will couple to the SM particles through higher dimensional operators, presumably 
mediated by particles of the dark sector which are somewhat heavier than the WIMP itself 
(and which may or may not carry SM gauge charges). 

This article is organized as follows. In Section [Til we write down a complete list of leading 
operators and analyze the constraints on the coefficients of these operators, assuming that 
only one operator is dominant at a time. In Section IHH we employ existing Tevatron and 
planned LHC searches to determine the constraints on the coefficients of the operators (or 
prospects for their discovery), respectively. Of the complete set of operators, some mediate 
substantial (i.e. not suppressed by the WIMP velocity) low energy WIMP-nucleus rates, and 
thus are constrained or may be discovered by direct detection experiments. These bounds 
and prospects are presented in Section HVl We conclude with comments on future directions 
in Section IVl 



II. EFFECTIVE FIELD THEORY OF WIMP COUPLINGS 

We consider the cases where the DM particle is a scalar or a fermion; if a scalar, it can 
be real or complex, and if a fermion, it can be Majorana or Dirac. Each of these cases 
is considered separately. We note that in principle, the WIMP could also be spin one or 
higher; we shall not consider these cases here since the couplings of such WIMPs are usually 
restricted by gauge invariance and other symmetries, and are more heavily model-dependent. 



3 



We shall be considering the situation where the WIMP (which we will generically denote 
X) is the only particle in addition to the standard model fields accessible to colliders. We 
will assume that x is °dd under some Z 2 symmetry (e.g. R-parity in supersymmetry, or KK- 
parity in extra dimensions), and hence each coupling involves an even number of WIMPs 
with the lowest dimensional operators we consider containing two WIMPs. We assume 
whatever particles mediate interactions between the WIMPs and the SM fields are somewhat 
heavier than the WIMPs themselves, with their leading effect manifest as higher dimensional 
operators in the effective field theory. For simplicity, we assume the WIMP is a singlet under 
the SM gauge groups, and thus possesses no tree-level couplings to the electroweak gauge 
bosons. We also neglect couplings with Higgs bosons. While the inclusion of such couplings 
in the effective theory is straightforward, we leave these cases for future work. Given the 
assumption that the WIMPs are SM singlets, the factor in each operator consisting of SM 
fields must also be invariant under SM gauge transformations. 

We note in passing that even for an electroweak singlet WIMP, the lowest dimensional 
operator linking, a pair of WIMPs to the SM fields contains two WIMPs and the SM Higgs 



bilinear |i?| 2 (27|. Such an interaction contributes to direct detection and collider processes 
involving WIMPs by inducing a x~X-h° interaction after electroweak symmetry-breaking 2 . 
While we do not consider this operator further, we note that for cases where the Higgs is 
heavy enough, it is effectively integrated out, leaving behind operators which we do consider. 

The next allowed class of operators have SM factors which are quark or lepton bilinears. 
The lepton bilinear couplings contribute only at a very suppressed level to direct detection 
or hadron collider production, leaving us with little to say about them. It would be very 
interesting to study constraints on such operators from indirect detection experiments and/or 
lepton colliders such as LEP-II. In this note we shall focus on quark bilinear operators of 
the form qTq, where T is a 4 x 4 matrix of the complete set, 

T = {l, 7 5 ,7 M ,7 M 7 5 ,^} ■ (1) 

We do not consider terms with derivatives acting on the quarks, which lead to higher di- 
mensional operators, more suppressed at low energies. 

Finally, we have operators mediating WIMP couplings to massless gauge fields. The 
leading operators are a magnetic moment coupling XG^xF^v and an electric dipole moment 



2 Collider and direct detection signals from this operator are explored in Ref. 28 1 
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coupling 'l&xFpv (which are only non- vanishing for a Dirac fermion WIMP), though 
given the unbroken U(1)em gauge invariance, they are likely induced at the loop level and 
thus may have small coefficients. Various ex per imental bounds and direct detection signals 



of these operators have been studied in Refs. 



29 



-38 



We do not consider collider constraints 



for these operators further here and leave astrophysics bounds for future work 39j . We also 
have couplings to G^Gap, where GG can either be a pair of electromagnetic or color field 
strengths, with gauge and Lorentz indices contracted in all possible ways to form a family 
of related operators. Here we focus on the operators involving color field strengths. Just as 
for quark operators, terms with derivatives acting on the gauge field strengths are higher 
order and more suppressed. 

All together, these higher- dimensional operators define an effective field theory of the 
interactions of singlet WIMPs with hadronic matter. It is expected to reasonably capture 
the physics provided the WIMP is somewhat lighter than the particles which mediate its 
interactions with the SM. It is a non-renormalizable field theory and thus must break down 
at some energy scale, represented by the masses of those particles which have been integrated 
out. The quantities M* which characterize the interaction strength of the interactions are 
functions of the masses and the coupling strengths of the mediating particles to WIMPs and 
SM fields, and can be computed in terms of the fundamental parameters for any specific UV 
theory of interest. 

What happens above the regime of validity of the effective theory depends on the UV 
completion, and is much more model dependent. Depending on the specifics of the UV 
theory, collider bounds may get stronger or weaker. For example, in supersymmetric theories 
our operators are UV completed into colored squarks which will be produced on-shell and 
may contribute to the jets + missing energy observable with larger rates than those we are 
computing in the effective theory. Other UV completions, such as a light neutral mediator, 
can lead to much weaker collider cross sections [9[, since far above the mediator mass the 
rate will fall with jet transverse energy as 1/-P t 2 , whereas in the effective theory the partonic 
reaction is flat with jet Pt, scaling as 1/Mjf. Thus, it should be borne in mind that our 
limits strictly speaking only apply when all mediator masses are much larger than the 
typical energy of the reaction, and in the absence of a picture of the UV theory, it is hard to 
know whether the bounds are over- or under-estimated when the effective theory description 
does not strictly apply. 
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For a given WIMP mass, there is a lower bound on M* such that one can imagine 
any weakly coupled UV completion. Since the operators mediate interactions with (at 
least) two colored SM fields coupled to two WIMPs, the simplest tree level UV completions 
have a single mediator particle and two interactions. The mapping to M* from the UV 
parameters thus involves an expression such as M* ~ Mj ^Jgxgh. where M is the mass of the 
exchanged particle, and g\ and g 2 are couplings. Since an effective theory description requires 
M > 2m x , and a preturbative theory gig 2 < (47r) 2 , a weakly coupled UV completion requires 
m x < 2ttM*, beyond which the UV completion becomes non-perturbative. In determining 
bounds, since there is no imaginable perturbative UV picture for m x < 27rM*, we cut off the 
bounded regions outside of this region of validity. Furthermore, for the effective theory to 
make sense, the mediator mass has to be larger than energy transfer through quarks at the 
collider environment. The limit, in which the effective theory breaks down, highly depends 
on the details of relevant patron energy and its distribution. Since M < 4nM* for the 
perturbative UV completion, our bounds are valid when the characteristic energy transfer 
is smaller than 4nM*. The detailed analysis of this limit is beyond the scope of this work, 
we will leave it for the future investigation. 

The coefficients of the operators are chosen to simplify comparisons to direct detec- 
tion experiments. For quark bilinears, the appropriate matrix elements (at low momentum 
transfer) are (N\m q qq\N) and (Nlq'j^qlN) which contribute to spin-independent scattering, 
(A r |g7 M 7 5 g| N), which contributes to spin-dependent scattering, and {N\qa^ u q\N), which 
couples to the magnetic moment of the nucleon. For the gluon operators, the relevant ma- 
trix element is (N\a s GG\N) . The scalar (and pseudo-scalar) quark bilinears are normalized 
by m q , which together with our choice of universal vector- type couplings has the added fea- 
ture of mitigating contributions to flavor changing processes from these operators, through 
the framework of minimal flavor violation 40|. For the gluon field strength operators, we 



normalize by a factor a s , which both anticipates their origin as loop processes and captures 
the dominant renormalization group evolution. The complete list of leading operators is 
given in Table I. The coefficients of these operators have been scaled by appropriate powers 
of M* (the value of which can be in principle different for each operator) to give the correct 
over-all dimension in the action. 



6 



Name 


Operator 


Coefficient 


Dl 


XXQQ 


m q /Ml 


D2 


X1 5 XQQ 


imjMl 


D3 


XXQ1 5 Q 




D4 


X1 5 XQ1 5 Q 


m q /Ml 


D5 


Xl^XQlfil 




D6 


X7 / Vx<Z7/Lt5 




D7 






D8 






D9 






D10 


xv^^xqvapq 


ijMl 


Dll 


xxG^G^ 




D12 


X^xG^ 




D13 


XXG^ 


ia s /±Ml 


D14 


Xl 5 xG„vG^ 


a s /4M3 



Name 


Operator 


Coefficient 


CI 


x ] xqq 


m q /Ml 


C2 


X f X<n 5 <? 




C3 


X f ^X97 M 9 


mi 


C4 


X t ^XQ7 At 7 5 9 




C5 


xhG^G^ 


a s /4M» 2 


C6 


X^G^G^ 


ia s /4M* 


Rl 


x 2 qq 


m q j1Ml 


R2 


x 2 qj 5 q 


irn q /2M* 


R3 


x 2 G^ 


a s /8M? 


R4 


x 2 G^ 


ia s /8Ml 



TABLE I: Operators coupling WIMPs to SM particles. The operator names beginning with D, C, 
R apply to WIMPS that are Dirac fermions, complex scalars or real scalars respectively. 

III. COLLIDER CONSTRAINTS 
A. Overview 

We can constrain M* for each operator in the table above by considering the pair pro- 
duction of WIMPs at a hadron collider: 



pp (pp) ->• XX + x. 



(2) 



Since the WIMPs escape undetected, this leads to events with missing transverse energy, 
recoiling against additional hadronic radiation present in the reaction. 

The most significant Standard Model backgrounds to this process are events where a Z 
boson decays into neutrinos, together with the associated production of jets. This back- 
ground is irreducible. There are also backgrounds from events where a particle is either 
missed or has a mismeasured energy. The most important of these comes from events pro- 



ducing W + jets, where the charged lepton from the W-dec&y is missed. Other backgrounds 
such as QCD multijet production (with the missing energy the result of mismeasuring the 
energy of one of more jets) are expected to be subdominant for the cuts chosen in the 



analyses 



41 



42] 



B. Tevatron Constraints 



The Tevatron has searched for signals of new physics with missing transverse energy in 
many channels. We will focus on monojet events, where the WIMPs recoil against a single 
jet, with restrictions on any additional SM radiation. We will compare the predictions of our 



431 ] . the events were required to 



effective theories with the results on monojet events from CDF |42|, |43j. We expect similar 
constraints can be derived from DO data, but choose to focus on the available CDF searches 
which utilize much greater integrated luminosity. In Ref. 
satisfy: 

• Events are required to have a leading Jet with transverse energy E t > 80 GeV; 

• Events must have net missing transverse energy E t > 80 GeV; 

• A second jet with E t < 30 GeV is allowed; 

• Events containing additional jets with E t > 20 GeV are vetoed. 



In order to simulate WIMP pair production events to 



compare to these bounds, we found 



44 



45] , where j is any parton other 



the partonic cross section for pp — > jxx using Comphep 
than the top quark, and is required to have E t > 80 GeV. At the parton level, this simulta- 
neously requires that Et > 80 GeV. We correct these parton-level estimates by an efficiency 
taking into account corrections from parton showering, hadronization, and energy smearing 
by the detector. This efficiency is computed by first hadronizing the generated parton- 
level events using Pythia 46|] (through the Comphep-Pythia interface i 
events are reconstructed at the detector level by passing them through PGS [48j] tuned to 
simulate the response of the CDF detector 3 and required to satisfy the detailed CDF anal- 
ysis cuts. The efficiency is defined as the ratio of the number of events after the PGS-level 



47l| ). The hadronized 



A previous study |6| found that this detector model was able to reproduce the backgrounds quoted in 



431 ] to the few percent level. 



cuts to the number at the parton-level. We find that efficiencies range from 30% to 48% 
for WIMPs with various spins and masses ranging from — 300 GeV. Given the relative 
insensitivity to the details, we choose a flat efficiency of 40% for all cases. 

We plot 2a lower limits on the scale of new physics M* for each operator as solid lines in 
figures I8~HTS| where for illustration, we also plot the lines resulting in the observed thermal 



relic density. Comparing with previous studies which analyzed D8 j^] and D 
and D8 (but with a different normalization between different flavors of quarks) 
rough agreement with these studies. It is worth noting that the CDF analysis was somewhat 
optimized for theories with large extra dimensions as opposed to xXJi an d it is possible that 
better bounds may be available for more optimized analysis strategies. 

C. LHC Prospects 

We also simulate the inclusive jets and missing transverse energy events at the LHC 



D4, D5, 
9], we find 



experiments for -y/i = 14 TeV and compare them with the analysis in Ref. 49j, which studied 
this signal in the context of a large extra dimensions search. In Ref. j^], the missing E t 
was required to be larger than 500 GeV and no veto on additional hadronic radiation was 
imposed. Additional hard isolated leptons were vetoed in order to reduce backgrounds from 
W + jets processes. Finally, the azimuthal angle between the the missing transverse energy 
and the second hardest jet was required to be A$ > 0.5, to mitigate the QCD background 
due to mismeasured jet energies. The expected number of background events after these 
cuts was simulated to be about B = 2 x 10 4 with 100fb _1 of data. 

Proceeding as before, we find the parton-level cross section <Jj XX (with E t of the jet greater 
than 500 GeV) using Comphep. We use Pythia and the PGS simulation with the generic 
LHC detector model to estimate an efficiency, defined as for the Tevatron study, of roughly 
80%. This is roughly in agreement with the efficiency for the mono jet signal from the large 
extra dimension model, which was found in Ref. 49| to be ~ 90%. 



We define the 5<r detection region for which the number of expected signal events 5* and 
background events B satisfy S/ \f~B > 5 for an integrated luminosity of 100 fb _1 . We plot this 
5cr reach as dashed lines in figures l8lfT6l Again, we note that this search was optimized for a 
large extra dimensions signal, and our knowledge of the LHC detector performance and SM 



backgrounds are expected to substantially improve over those available to Ref. 49j. Given 
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our positive results, a dedicated reanalysis by the collaborations would be very interesting. 



IV. IMPLICATIONS FOR DIRECT DETECTION 

Our new bounds on the strength of interactions of WIMPs with hadrons can be translated 
into constraints on the possible contributions to direct detection cross sections for each of 
those interactions. Only some operators contribute to direct detection in the limit of low 
momentum transfer, and the remaining operators are suppressed by powers of the WIMP 
velocity, generically expected to be of order ~ 10~ 3 . For each contributing operator we 
employ the expectation value of the partonic operator in the nucleon 50||. This, combined 
with the kinematics of WIMP-nucleon scattering, results in cross sections 



Dl 



cr 



(J, 



D5,C3 




1.60 x 10 



- 37 cm 2 



Jbc_\ 
lGeV 



2 /20GeVy 



1.38 x l(T 37 — 2 



cm 



the ) 
lGeV 



D8,D9 







9.18 x 10" 40 cm 



i 2 (J*J\ 
VlGeV 



n 

2 / 300GeV \ 4 
) 1 M* ) 



300GeV\ 4 



cr, 



Dll 



3.83 x 10" 41 cm 



2 ( \ 
VlGeV 



/m 8 



_C1,R1 




2.56 x 10" 36 cm 2 



H X \ 2 f 10GeV\ 2 / lOGeVV 
lGevJ V m x ) V M * ) 



cr. 



C5,R3 







7.40 x 10" 39 cm 



( the \ 
VlGeV 



2 /10GeV\ 2 /60GeV\' 



) \ m x J { M m J 



(3) 
(4) 
(5) 
(6) 
(7) 
(8) 



where [i x is the reduced mass of the WIMP-nucleon system. 

The behavior at low WIMP masses is affected strongly by the spin of the WIMP itself. 
For a fermion WIMP the direct detection cross section for a fixed coupling is largely flat 
as the WIMP mass is decreased, until the WIMP is lighter than the nucleon mass. For 
a scalar WIMP, except for the vector-type interaction C3, the mass appears explicitly in 
the expression for the cross section, causing the cross section for smaller WIMP masses to 
notably increase (provided m x < M*). This has the net effect of weakening the impact of 
the collider bounds on the direct detection parameter space for of very light scalar dark 
matter with respect to those for fermion dark matter. 

We also notice that the collider bounds on direct detection of non-self-conjugate fields 
are stronger by a factor of 2 in cross section than those on self-conjugate fields. This is an 
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expected result, as the phase space for direct detection is unchanged by this factor but the 
phase space of collider production is suppressed by a factor of 2 for self-conjugate fields. 

We plot the effective collider constraints on the WIMP direct detection parameter space 
in figures \T§7\ including the most relevant direct detection constraints for comparison. We 
see that in all cases, colliders can probe regions of very light WIMP masses more effectively 
than direct detection experiments, which become limited by energy thresholds for extremely 
light WIMPs. Indeed, for many operators, the direct detection rates are expected to be very 
small because of the velocity suppression, and colliders become the only way to effectively 
probe WIMP-hadron interactions. In the case of a WIMP whose dominant recoil is through 
a spin- dependent interaction, collider constraints are already much stronger even than the 
expected reaches of near-future direct detection experiments. Thus, if such an experiment 
were to observe a positive signal, the collider constraints would immediately imply a break- 
down of the effective field theory at collider energies, revealing the existence of a light 
mediator particle. 

V. CONCLUSIONS AND OUTLOOK 

We have studied constraints on dark matter models coming from collider experiments, 
specifically from extant searches at the Tevatron and future searches at the LHC. This 
extends our previous work where we considered the case of dark matter particles which were 
Majorana fermions to the cases where the dark matter is a Dirac fermion or a scalar (either 
real or complex). 

Our results are qualitatively similar to our previous paper. In general collider constraints 
are very strong for lighter dark matter and fall off when the dark matter mass exceeds the 
typical energy reach of the collider. The constraints also depend on the coupling of the dark 
matter; if the dark matter primarily couples to gluons, the constraints from colliders become 
especially strong. 

One of the most interesting results is that collider constraints on spin dependent interac- 
tions are stronger than direct searches over a significant portion of parameter space. In the 
event that direct searches find a signal in this region while no signal is found at colliders, 
this will suggest that dark matter is a WIMP of spin 1 or higher, or that there exists a light 
mediator particle UV completing the interaction operators in such a way as to weaken the 
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?IG. 1: Current experimental limits on spin- independent WIMP direct detection from CRESST 
53], Xenon 10 Q, CoGeNT[3], and Xenon 100 Q, (solid lines ja labeled), as well 
JeNT favored region Q 



5J, CDMS 



551 ] and Xenon 100 



and future reach estimates for SCDMS 
5a ]. where we have chosen the line using a threshold of 3PE and the conservative extrapolation 
of £ e ff (dashed lines as labeled). Also shown are the current Tevatron exclusion for the operator 
Dll (solid magenta line) as well as LHC discovery reaches (dashed lines as labeled) for relevant 
operators. 



collider bounds. The case of a light mediator with a particular 

dark matter + dark matter -h- SM-neutral mediator -h- SM + SM 

completion structure was considered in {9]. Beyond these particular constructions, many 
models have additional light states which UV complete the interactions between the dark 
matter and the Standard Model through a 

dark matter + SM -H- SM-charged mediator -H- dark matter + SM 

topology. It would be relatively simple to consider a complete set (as dictated by SM gauge 
and Lorentz invariance) of UV completions, and it would be interesting to see how our 
bounds are modified in the presence of such new states, and whether new collider signals 



12 




10"' 



10 



m, (GeV) 



10 z 



KIMS p limits 



10 3 



FIG. 2: Current experimental limits on spin-dependent WIMP direct detection from Picasso 



KIMS 



58|, and Xenon 10 



591 ] . Also shown are the 



54( | , as well as the future reach of DMTPC 
current Tevatron exclusions (solid lines as labeled) and LHC discovery reaches (dashed lines as 
labeled) for relevant operators. 



can be found to constrain such models. We leave detailed exploration of these issues for 
future work. 

Finally, we note that while effective theories may not always capture our favorite parame- 
ters of our favorite UV-complete models, they do provide a language to describe WIMP-SM 
interactions which captures a wide class of theories in a fairly model-independent fashion. 
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FIG. 3: Same as Fig.[H but with Tevatron exclusions and LHC reaches for complex scalar WIMP 
operators, as labeled. 




FIG. 4: Same as Fig.[H but with Tevatron exclusions and LHC reaches for real scalar WIMP 
operators, as labeled. 
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FIG. 5: A cartoon representation of previous limits due to direct detection experiments as well as 
constraints from Earth heating and cosmic rays 5l[, with new exclusions from Tevatron searches 



for Dirac WIMPs superimposed. 




FIG. 6: Same as Fig.[5l but with complex scalar WIMP limits superimposed. 
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FIG. 7: Same as Fig. [SI but with real scalar WIMP limits superimposed. 
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FIG. 8: Constraints on the operators Dl-4. Solid lines are Tevatron 2a lower limits, dashed lines 
are LHC 5a discovery reach lower limits, and dot-dashed lines indicate the values of M* necessary 
for the WIMP to have the correct relic abundance in absence of any other interactions. The 
curves for operators Dl and D2 are largely degenerate with those for D3 and D4, respectively. 
The gray filled region indicates where the effective field theory breaks down, possessing no simple 
perturbative UV completion. 
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FIG. 9: Same as Fig. [51 but for the operators D5 and D6 which are largely degenerate with D7 and 
D8, respectively. 
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FIG. 10: Same as Fig.[8j but for the operators D9 and D10. 
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FIG. 11: Same as Fig. El but for the operators Dll and D12 which are largely degenerate with D13 
and D14, respectively. 
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FIG. 13: Same as Fig. [51 but for the largely degenerate operators C3 and C4. 
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FIG. 16: Same as Fig.0 but for the largely degenerate operators R3 and R4. 
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